Teleportation fidelities of squeezed states from thermodynamical state space measures 
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We introduce a "microcanonical" measure (complying with the 'general canonical principle') over the sec- 
ond moments of pure bosonic Gaussian states under an energy constraint. We determine the average fidelity 
for the teleportation of states distributed according to such a measure and compare it to a threshold obtained 
from a feasible classical strategy. Furthermore, we show that, under the proposed measure, the distribution of 
the entanglement concentrates around a finite value at the thermodynamical limit and, in general, the typical 
entanglement of Gaussian states with maximal energy E is not close to the maximum allowed by E. 
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Introduction - Besides having been at the core of theoretical 
and experimental quantum optics right from its early stages, 
Gaussian states have recently acquired a major role in quan- 
tum information science, in the so called 'continuous variable' 
(CV) scenario |Q]] . Indeed, some of the most spectacular im- 
plementations of quantum information protocols to date are 
based on Gaussian states, with the prominent example of de- 
terministic teleportation yfl. In the analysis of most such im- 
plementations, the proper assessment of figures of merit re- 
quires the average over a distribution (a "measure") of states 
in which input quantum information is encoded. For instance, 
in the case of quantum teleportation of coherent states, the 
theoretical average fidelity (between input and output states) 
is determined by assuming a particular distribution of input 
coherent states yfl. In the present paper we propose a mea- 
sure on the set of pure Gaussian states, whose introduction 
will be thoroughly motivated by fundamental statistical ar- 
guments 0]. We shall focus on the second moments of the 
quadrature operators (while the measure usually employed to 
analyze teleportation of coherent states ^ essentially encom- 
passes first moments), covering the whole set of pure Gaus- 
sian states with null first moments. As we will mention later, 
first moments may be accomodated as additional variables in 
the presented framework. 

The importance of determining a suitable measure over a 
set of states is not merely a theoretical issue, as the evaluation 
of classical thresholds for the figures of merit is crucial in es- 
tablishing whether practical realizations of quantum protocols 
actually out-perform competing classical strategies J^. We 
shall thus apply the proposed measure to determine the aver- 
age teleportation fidelity of pure Gaussian states with varying 
second moments, and shall compare such a fidelity to a cor- 
responding "classical" threshold. Moreover, to further illus- 
trate the potentialities of a measure on second moments, we 
will address the "typical" entanglement J5[] of pure Gaussian 
states under an energy constraint. The very construction of 
the measure will imply that the distribution of the von Neu- 
mann entropy of any finite subsystem 'concentrates', both at 
the thermodynamical limit and for finite numbers of modes, 
around a finite 'thermal' average, well away from the allowed 
maximum. 

Preliminary facts and notation - We consider continuous 



variable quantum mechanical systems described by n pairs 
of canonically conjugated operators {xj, Pj} with continuous 
spectra. Grouping the canonical operators together in the vec- 
tor R = (xi, . . . ,x n ,pi, . . . ,p n ) T allows one to express the 
canonical commutation relations as [Rj, R^) = 2i Slj k , where 
the 'symplectic form' SI has entries Sljk = (Jj+n.fc — $j,k+n 
for j, k = 1, . . . ,2n. Any state of an n-mode CV sys- 
tem is described by a positive, trace-class operator g. For 
any state g, let us define the 2n x 2n matrix of second mo- 
ments, or "covariance matrix" (CM), <x with entries a^h 
Tr [{Rj,R k }Q]/2 - Tr [RjQ]Ti [R k g], In the following, we 
will refer to the 'energy' of a state g as to the expectation value 
of the operator Hq = 2j=i {%j + Pj) ( n °te that, in our con- 
vention, the vacuum of a single mode has energy 2). This defi- 
nition corresponds to the energy of a free electromagnetic field 
in the optical scenario (and to decoupled oscillators in the gen- 
eral case). Neglecting first moments, the energy is determined 
by the second moments according to Tr (gHo) = Tr (er). 

Gaussian states are defined as the states with Gaussian char- 
acteristic functions and quasi-probability distributions. All 
pure Gaussian states can be obtained by transforming the vac- 
uum under unitary operations generated by polynomials of the 
second order in the canonical operators. Operations gener- 
ated by first order polynomials in the quadratures correspond 
to local displacements in the first moments, and will thus be 
disregarded. As for second order transformations, they can 
be mapped into the group Sp2n.n of real symplectic transfor- 
mations, by virtue of the so called metaplectic representation 
(recall that S G SL(2n, R) : S G Sp 2n ,& ^ S T SLS = SI). 
As a consequence of such a mapping, the CM er of any pure 
Gaussian state can be written as er = S T S JUS- 
General canonical principle and microcanonical measure. - 
We will now proceed to define a measure over the set of pure 
Gaussian states, which will be referred to as 'microcanoni- 
cal' (for reasons which will be clear shortly). Henceforth, the 
shorthand notation x will stand for the average of the quan- 
tity x with respect to such a measure. Since we will adopt a 
constructive approach, based on the gradual enforcement of 
specific conditions on the measure, the notation x will appear, 
with no ambiguity, before the definition of the measure itself. 

Because the symplectic group is non-compact, an invariant 
Haar measure on the whole group (from which a measure for 
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the second moments of pure Gaussian states could be derived 
via the equation er = S 1 S) would be non-normalizable (and 
thus "unphysical", giving rise to distributions with unbounded 
statistical moments). The first natural prescription to tame the 
non-compact nature of the group consists in the introduction 
of a constraint on the total energy of the system, which we 
will denote by E (hence the designation "microcanonical" at- 
tached to the measure). Even so, no natural invariant measure 
emerges. However, let us recall that an arbitrary symplectic 
transformation S can be decomposed as S = O' (Z © Z )0, 
where O, O' G K(n) = Sp(2n,H) D SO(2n) are orthogo- 
nal symplectic transformations, while Z is an n x n diagonal 
matrix with eigenvalues Zj > 1 for 1 < j < n $7]. The 
set of transformations of the form Z © Z^ 1 is a non-compact 
subgroup of Sp2n.~R (corresponding to local squeezings). The 
virtue of such a decomposition, known as the 'Euler' decom- 
position, is immediately apparent, as it allows one to distin- 
guish between the degrees of freedom of the compact sub- 
group ('angles', collectively denoted by $, which do not affect 
the energy) and the degrees of freedom Zj's with non-compact 
domain. In particular, applying the Euler decomposition to the 
CM er of generic pure states leads to er = O 1 {Z 2 © Z~ 2 )0. 
Moreover, we recall that the compact subgroup K(n) is iso- 
morphic to U(n) J3]. As dictated by the Euler decomposi- 
tion, we assume the n 2 parametres & of the transformation 
O G K(n) to be distributed according to the Haar measure of 
the compact subgroup K{n), which can be carried over from 
U(n) through the isomorphism recalled above and whose in- 
finitesimal element will be denoted by dp,H ($)■ 

We are thus left with the parameters Zj's alone, for which 
a 'natural' measure has not yet emerged. To constrain the 
choice of such a measure, we will invoke a fundamental statis- 
tical argument. In their alternative, 'kinematical' approach to 
statistical mechanics, Popescu et al. (01 define a general prin- 
ciple, which they refer to as general canonical principle, stat- 
ing that "Given a sufficiently small subsystem of the universe, 
almost every pure state of the universe is such that the subsys- 
tem is approximately in the 'canonical state' q c ." The 'canon- 
ical' state q c is, in our case, a Gaussian thermal state, with 
CM er c = (1 + T/2)t. Here the 'temperature' T is defined 
by passage to the thermodynamical limit, that is for n — > oo 
and E — ► oo, (E — 2n)/n — ► T (assuming ks = 1 for the 
Boltzmann constant). For ease of notation, in the following, 
the symbol ~ will imply that the equality holds at the ther- 
modynamical limit, e.g.: (E — 2n)/n ~ T. Because the state 
q c is Gaussian with null first moments, the general canonical 
principle can be fully incorporated into our restricted (Gaus- 
sian) setting. To this aim, let us recast the principle in terms of 
mathematical conditions to be fulfilled by the underlying mea- 
sure on pure Gaussian states. Recall that partial tracing (ob- 
viously a Gaussian operation) amounts, at the level of CM's, 
to simply pinching the submatrix of er pertaining to the rele- 
vant modes. Then, in order for the measure to comply with 
the general canonical principle, one has to require 



The second condition rephrases the prescription "almost ev- 
ery pure state", requiring 'concentration of measure' at the 
thermodynamical limit (in fact, in conjunction with the first 
equation, it implies that the variance of the entries of the CM 
vanishes at the thermodynamical limit). The previous condi- 
tions, which are highly desirable to single out a measure nat- 
urally endowed with physical and statistical significance, will 
greatly restrict the possible choices for the distribution of the 
variables Zj's. 

In order to show this, we first work out the averages of the 
entries of er over the Haar measure of the compact subgroup. 
This task can be accomplished relying only on some basic 
properties of the integration over the unitary group, derived 
from simple symmetry arguments (see lioll . a detailed deriva- 
tion can be found in Jg]), leading to 
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The convenience of a parametrisation through the variables 
Ej = (z 2 + zj 2 ), representing the local energies of the decou- 
pled modes, is now apparent. The same arguments, based on 
symmetry and normalization (reflecting orthogonality), can be 
applied to the average a 2 k , leading to 
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(l + T/2) 2 S jk . (1) 



Now, the desired agreement of Eqs. ( 1213b with Eq. ([T]i single 
out a restricted class of measures for the variables Ej's. Most 
notably, any measure such that the local energies {Ej} are, 
at the thermodynamical limit, independent, identically dis- 
tributed ("i.i.d.") variables with average E = (T + 2) com- 
plies with the general canonical principle |Q] . In point of fact, 
at the thermodynamical limit, only the averages EjE k with 
j 7^ k matter in the computation of the variance, as their num- 
ber scales as n 2 , while the number of terms in E 2 is clearly 
linear in n, and their contribution gets suppressed by the fac- 
tor 1/n 2 (deriving from the integration over the Haar measure 
of a term of degree two in the compact transformations' en- 
tries). The same argument holds for the square of the quantity 

2 

WJk of Eq. (f2]). For i.i.d. variables, Ej ~ E and EjE k — E 
Vj 7^ k, thus implying the vanishing of the variance at the 
thermodynamical limit. 

To complete the definition of the measure, we have to spec- 
ify a distribution of the Ej "s in agreement with the previous 
requirements. Recovering the energy constraint E, we will as- 
sume a Lebesgue ('flat') measure for the local energies Ej's 
inside the region r# = {E : \E\ < E}, bounded by the lin- 
ear hypersurface of total energy E (here, E = (Ei, . . . , E n ) 
denotes the vector of energies, while \E\ = Y^j=i Ej). More 
explicitly, denoting by dp(E) the probability of the occur- 
rence of the energies E and by V = (E — 2n) n /n\ the 
volume of the region Te, one has dp(E) = d n E/V = 
(dE x .. .dE n )/V if E G T E and dp(E) = otherwise. 
Notice that such a flat distribution is the one maximising the 
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entropy in the knowledge of the local energies of the decou- 
pled modes. In this specific sense such variables have been 
privileged, on the basis of both mathematical (Euler decom- 
position and Haar averaging over the compact subgroup) and 
physical (general canonical principle and analogy with the mi- 
crocanonical ensemble) grounds. In full analogy with the 
equivalence of statistical ensembles, the Ej's become i.i.d. 
at the thermodynamical limit. In fact, the marginal density 
of probability P n (Ej , E) for each of the energies Ej given 

by P n (Ej,E) = (l - |^|)" . At the thermody- 

namical limit, the upper integration extremum for each Ej 

Ej-2 

diverges and P n (Ej,E) — > c ~ /T, so that the decou- 
pled energies are distributed according to independent Boltz- 
mann distributions with the parameter T playing the role of 
a temperature, and their averages satisfy Eq. ([T). The micro- 
canonical measure is thus consistent with the general canon- 
ical principle |@]. The 'microcanonical' average Q over pure 
Gaussian states at energy E of the quantity Q(E,d) deter- 
mined by the second moments alone will thus be defined as 
Q = Af J du ff ($) J T dEQ(E, where the integration 
over the Haar measure is understood to be carried out over 
the whole compact domain of the variables fl. 
Typical entanglement - Here, we concisely address the sta- 
tistical properties of the bipartite entanglement of pure Gaus- 
sian states under the microcanonical measure [6]. Let us con- 
sider the von Neumann entropy S of the reduced state g m 
of a finite number of modes m with CM 7, thus quantify- 
ing its entanglement with the remaining (n — m) modes of 
the globally pure state. At the thermodynamical limit the dis- 
tribution of the CM 7 concentrates, with vanishing variance, 
around a thermal state with CM (T/2 +1)1, according to 
the general canonical principle. Therefore, the distribution 
of the von Neumann entropy of the reduction S concentrates 
around the von Neumann entropy of a thermal state jfj. In 
formulae: S ~ m/(l + T/2) and (]P - S 2 ) ~ 0, where 
/(.*) ee (x+l)log 2 [(x+l)/2]/2-(x-l)log 2 [(x-l)/2]/2. 
Notice that the maximal local von Neumann entropy of any 
(finite or infinite) subsystem diverges at the thermodynamical 
limit (as, in principle, all the infinite energy could be concen- 
trated in only two modes - owned by the two distinct sub- 
systems -, thus yielding an infinite entropy for each subsys- 
tem). For finite n, the microcanonical measure is apt to be 
investigated numerically by direct sampling, allowing one to 
study the distribution of the von Neumann entropy for differ- 
ent m, n and E. Even for small n - well before the onset of 
thermodynamical concentration of measure around the finite 
thermal average - the entanglement of pure Gaussian states 
distributes, for small enough energies, around values gener- 
ally distant from the finite allowed maximum (e.g., for m = 1 
and E = lOn, the difference between the maximum and the 
average S is, respectively, 4.0 and 13.6 standard deviations for 
n = 5 and n = 20). This is at striking variance with results 
obtained, adopting different measures, in finite dimensional 
systems 0]. The equipartition of energy, imposed by the gen- 
eral canonical principle, prevents the entanglement of finite 
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FIG. 1: Average teleportation fidelity for r = 1 (blue curve) and 
r — 0.5 (green curve) versus the heterodyne threshold (red dots), as 
functions of the maximal energy E. 



subsystems from concentrating around the maximum. Notice 
that the concentration of measure would also occur for a distri- 
bution of states with fixed total energy E (and all the variables 
{Ej } still Lebesgue-distributed under such a constraint). This 
is the case as such a measure is equivalent, at the thermody- 
namical limit, to the one we are considering (they both con- 
verge to exactly the same Boltzmann distribution). 
Teleportation fidelities - Let us now consider a practical situ- 
ation, in which two parties (Alice and Bob) want to commu- 
nicate through quantum teleportation. Instead of contenting 
themselves with coherent states, Alice and Bob are interested 
in exchanging single-mode pure Gaussian states with arbitrary 
CM (wherein the quantum information is encoded). To this 
aim, they employ the usual CV teleportation scheme, based 
on homodyne measurements and on the sharing of a two-mode 
squeezed state with squeezing parameter r (essentially quan- 
tifying the entanglement exploited in the teleportation pro- 
cess, see for a detailed description of the scheme). Sup- 
pose, quite reasonably, that Alice generates (and sends to Bob) 
states with a flatly distributed energy up to a maximal value 
E and random optical phase, whose distribution is thus de- 
scribed by the microcanonical measure. The microcanonical 
average fidelity T (defined as the average, over the distribu- 
tion of input states and of measurement outcomes, of the over- 
lap I (in I out) 1 2 between the input state \in) and output state 
\out)) can be straightforwardly determined II ill , and found to 
2e 2r (Vl + c- 
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To properly assess the effectiveness of the standard teleporta- 
tion protocol in transmitting second moments, let us compare 
the previously obtained fidelity to an appropriate "classical 
threshold" T c i (as customary in the literature on teleportation, 
"classical" refers to a procedure where no entanglement is ex- 
ploited). The kind of classical strategy we will consider is 
as follows: Alice measures her mode by heterodyne detec- 
tion (corresponding to the positive operator valued measure- 
ment with elements |a)(a|/7r, where \a) is a coherent state 
1 12J) and sends her result (a complex number a, represent- 
ing the heterodyne signal) to Bob who reproduces a centered, 
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pure state (belonging to the original distribution) lfl3ll . Bob's 
choice, depending on a, has been optimised numerically and 
the resulting fidelity has been averaged over the input distri- 
bution of states, to obtain the "heterodyne" classical threshold 
T c i (lower bounding the actual classical threshold). For suffi- 
ciently small E, one can also approximate such a threshold - 
obtained numerically - as T c i ~ 1 — fc arcsinh(\/.E — 2), for 
k = 0.317576 (such a fit is reliable within 0.002 in the range 
2 < E < 8). Comparing the experimental average fidelity to 
the previous formula would tell Alice and Bob whether their 
precious quantum entanglement is offering an actual advan- 
tage over a viable, 'cheaper' protocol based on disjoint mea- 
surements and reconstructions of the states 11411 . Fig. [TJ shows 
that, for a given r, the classical strategy beats CV teleporta- 
tion for small enough E. Actually, this fact is not surprising: 
it simply results from the inadequacy of the standard telepor- 
tation protocols when the input alphabet is overly restricted, 
and occurs in the teleportation of coherent states as well, if 
the choices of the coherent amplitudes are sufficiently con- 
strained (the reader might think about the limiting instance for 
which the vacuum is the only input state: then the teleporta- 
tion protocol, completely based on probabilistic measurement 
outcomes, fails to yield a fidelity equal to one, whereas the 
classical protocol is set to always return the vacuum in such 
an instance). However, there always exists a value of E above 
which the C V protocol starts outperforming the classical strat- 
egy. For instance, in the experimentally realistic case r « 1, 
such a threshold is remarkably low, being around E w 2.16. 
On the other hand, for any value of E, the teleportation proto- 
col may always exceed the classical threshold for high enough 
r. Clearly, for any finite E, one has lim^oo T = 1 (tele- 
portation with unlimited resources is always perfect). Also, 
numerics unambiguously show that lims^oc T c i = 0. Like- 
wise however, for any finite r, one has \\mE-*oo T — 0. This 
limiting behaviour is quite remarkable and might inspire fu- 
ture inspection into the matter: when the alphabet of states is 
enlarged to encompass all the possible second moments, the 
fidelity of the standard teleportation protocol vanishes (as op- 
posed to what happens for coherent states, where the fidelity 
stays constant even if the alphabet is extended over an un- 
bounded domain). This suggests that, possibly, a modified 
protocol where Bob can act unitarily on the second moments 
could grant better fidelities when the teleportation of second 
moments is concerned. 

Outlook - The study of generic entanglement and of figures 
of merit for teleportation are only examples of the potential 
applications of the microcanonical measure. For instance, 
the compliance with the general canonical principle renders 
the measure suitable to describe the thermaHzation of systems 
in dynamical situations lfl5ll . Relating the measure to distri- 
butions derived from a randomizing process (in the spirit of 
Ref. IU6I1 ) is a further line of development opened up by the 
present investigation. 
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